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I am repeating the proposal abstract here to provide the context for 
the final report: 

X propose to continue providing observers with basic data for interpreting 
spectra from stars, novas, supernovas, clusters, and galaxies. These data 
will include allowed and forbid den line lists, both laboratory and computed, 
for the first five to ten ions of all atoms and for all relevant diatomic 
molecules. X will eventually expand to all ions of the first thirty elements 
to treat far UV and X-ray spectra, and for envelope opacities. X also include 
triatomic molecules provided by other researchers. X have made CDs with 
Partridge and Schwenke's water data for work an X stars. The line data also 
serve as input to my modal atmosphere and synthesis programs that generate 
energy distributions, photometry, limb darkening, and spectra that can be used 
for planning observations and for fitting observed spectra. The spectrum 
synthesis programs produce detailed plots with the lines identified. Grids 
of stellar spectra can be used for radial velocity-, rotation-, or abundance 
templates and for population synthesis. X am fitting spectra of bright stars 
to test the data and to produce atlases to guide observers. For each star 
the whole spectrum is computed from the UV to the far XR. The line data, 
opacities, models, spectra, and programs are freely distributed on CDs and on 
my Web site and represent a unique resource for many N ASX programs. 

X am now in full production of new line lists for atoms [, except that 
X have not been getting paid since May 15, 2004, and except that X am going 
to run out of money for supplies and equipment] • X am computing all ions of 
all elements from H to Zn and the first 5 ions of all the heavier elements, 
about 800 ions. For each ion X treat as many as 61 even and 61 odd 
configurations, computing all energy levels and eigenvectors • The 
Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least 
squares fitting the observed energy levels. The average energy of each 
configuration is used in computing scaled -Thomas -Fermi -Dirac wavef unctions 
for each configuration which in turn are used to compute allowed and forbidden 
transition integrals. These are multiplied into the LS allowed and forbidden 
transition arrays. The transition arrays are transformed to the observed 
coupling to yield the allowed and forbidden line lists. Results are put on 
the web as they are finished. [Provided X get funding,] there will be more 
than 500 million lines. X will then compare ion by ion, to all the laboratory 
and computed data in the literature and make up a working line list for 
spectrum synthesis and opacity calculations with the best available data. 

As the laboratory spectrum analyses are improved, X will redo the calculations 
with the new energy levels. 

My original plan whan X started the new calculations was to run through 
all the atoms using my old Cray programs from the 1980's that were limited to 
1100 x 1100 arrays in the Hamiltonian for each J. Then X would go back and 
rerun the more complicated cases with 3000 x 3000 arrays so that X could 
include many more configurations and many more configuration interactions . 

At present X am limited to 61 even and 61 odd configurations and X try to 
include everything up through n * 9. The current program runs on Alpha 
workstations. X decided to test the big program on Fe I and Fe XX to see 
whether there was any great difference in the low configurations compared to 
those from the Cray program. 
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Besides increasing the number of El lines by a factor of 6 to 7.7 
million, there was an unexpected result: the electric quadrupole transitions 
were 10 times stronger than before because the transition integrals are 

weighted by r A 2 they become very large for high n, and because there are 

numerous configuration interactions that mix the low and high configurations. 
As a check X was able to reproduce Garstang's (1962) lower results by running 
his three configurations with my program. Since my model atom is still only 
a subset of a real Fe IX ion, the true quadrupole A values are probably 
larger than mine. The magnetic dipole lines are affected by the mixing but 
the overall scale does not change. 

Because of this discovery X decided that there was no point in computing 
the small array cases. X have been running with as many configurations as 
X can and with thousands of parameters in the Hamiltonian. The computer 
runs take much longer to set up and produce than X had expected. X have 
concentrated on redoing the low iron group spectra, especially to get data 
for supernova modelers. X have done only Ca X — Zn X, Ca XX — Zn II, 

Cu I — Cu XXXX, Zn X - Zn XXX, for practice at high stages of ionization, 

C X, C XX, S X, and Cl X and Ag X for people who were working on the 
laboratory spectra. Check my web site kurucz.harvard.edu for current 
additions . 

My latest calculations have been for carbon X and sulphur X, and 
silicon X is under way using the same elaborate approach as for C I, 
which took many months to do. These line lists greatly increase the 
number of lines in the ultraviolet, in the visible, and especially in 
the infrared. They will increase the opacity in A, F, and G stars. 

They will account for many unidentified lines in the sun. Here are 
parts of the .README files: 

0600 .README 9 October 2004 0600 C X 

Carbon X has now been completely redone and expanded and put in 
directory /atoms/0600. The computation was divided into four sections: 
suffix X, roughly all levels up through n * 9-12 with all configuration 
interactions; suffix R for Ritz, all higher levels up though n « 20 for 
low 1 with no configuration interaction; suffix P for polarization, all 
higher levels up though n * 20 for high 1 with no configuration 
interaction; and, by extrapolation, levels up through n » 99 that 
connect to s2p2, sp3, 2p3s, 2p3p, 2p3d levels and merge into a continuum. 

All higher levels up to n ■ 20 were filled in by Ritz expansion or the 
polarization formula for average energies. The polarization average 
energies are also given for all configurations with 1<21 for all n<100. 

These are the configurations in each section. 

X 


1 2p2 

2 2p3p 

3 2p4p 

4 2p5p 

5 

2p6p 

6 2p7p 

7 2p8p 

8 2p9p 

9 2pl0p 

A 2pllp 

B 

2p4 

C 2p3 3p 

D 2p3 4p 

E 2p3 5p 

F 2p3 6p 

G 2p3 7p 

H 

2p3 8p 

X 2p3 9p 

J 2p3 10p 

K 2p3 lip 

L 2p4f 

M 2p5f 

N 

2p6f 

O 2p7f 

P 2p8f 

Q 2p9f 

R 2pl0f 

S 2pllf 

T 

2p3 4f 

U 2p3 5f 

V 2p3 6f 

W 2p3 7f 

X 2p3 8f 

Y 2p3 9f 

Z 

2p3 lOf 

a 2p3 Ilf 

b 2p6h 

c 2p7h 

d 2p8h 

e 2p9h 

f 

2pl0h 

g 2p8k 

h 2p9k 

i sp2 3s 

j sp2 4s 

k sp2 5s 

1 

sp2 6s 

m sp2 7s 

n sp2 8s 

o sp2 9s 

p sp2 108 

q sp2 11s 

r 

sp2 3d 

s sp2 4d 

t sp2 5d 

u sp2 6d 

v sp2 7d 

w sp2 8d 

X 

sp2 9d 

y sp2 lOd 


z sp2 lid 
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1 2p3s 

2 2p4s 

3 2p5s 

4 2p6s 

5 2p7s 

6 

2p8s 

7 2p9s 

8 2pl0s 

9 2plls 

A 2pl2s 

B 2p3 3s 

C 

2p3 4s 

D 2p3 5s 

E 2p3 68 

F 2p3 7s 

G 2p3 8s 

B 2p3 9s 

I 

2p3 108 

J 2p3 11s 

K 2p3 12s 

L 2p3d 

M 2p4d 

N 2p5d 

o 

2p6d 

P 2p7d 

Q 2p8d 

R 2p9d 

S 2pl0d 

T 2plld 

u 

2pl2d 

V 2p3 3d 

W 2p3 4d 

Z 2p3 5d 

T 2p3 6d 

Z 2p3 7d 

a 

2p3 8d 

b 2p3 9d 

c 2p3 lOd 

d 2p3 lid 

• 2p3 12d 

£ 2p5g 

g 

2p6g 

b 2p7g 

i 2p8g 

j 2p9ff 

k 2pl0g 

1 2p7i 

a 

2p8i 

n 2p9i 

o 2p91 

p sp3 

G sp2 3p 

r sp2 4p 

8 

sp2 5p 

t sp2 6p 

u 8p2 7p 

▼ sp2 8p 

w sp2 9p 

3C sp2 lOp 

y 

sp2 lip 

z sp2 12p 







R 

1 2pl2p 

2 2pl3p 

3 2pl4p 

4 2pl5p 

5 2pl6p 

6 

2pl7p 

7 2pl8p 

8 2pl9p 

9 2p20p 

A 2pl2f 

B 2pl3£ 

C 

2pl4£ 

D 2pl5£ 

E 2pl6f 

F 2pl7£ 

G 2pl8£ 

B 2pl9£ 

I 

2p20£ 

J 2pllh 

K 2pl2h 

L 2pl3b 

M 2pl4b 

H 2pl5h 

O 

2pl6b 

P 2pl7h 

Q 2pl8h 

R 2pl9b 

S 2p20b 

T 2pl0k 

u 

2pllk 

V 2pl2k 

W 2pl3k 

Z 2pl4k 

T 2pl5k 

Z 2pl6k 

a 

2pl7k 

b 2pl8k 

c 2pl9k 

d 2p20k 

e 

£ 

g 


b 

i 

i 

k 

1 

a 


n 

o 

P 

q 

r 

8 


t 

u 

▼ 

w 

X 

y 


Z 

1 2pl3s 

2 2pl4s 

3 2pl5s 

4 2pl6z 

5 2pl7s 

6 

2pl8s 

7 2pl9z 

8 2p20s 

9 2pl3d 

A 2pl4d 

B 2pl5d 

C 

2pl6d 

D 2pl7d 

E 2pl8d 

F 2pl9d 

G 2p20d 

B 2pllg 

X 

2pl2g 

J 2pl3g 

K 2pl4g 

X. 2pl5g 

X 2pl6g 

V 2pl7g 

o 

2pl8g 

P 2pl9g 

Q 2p20g 

R 2pl0i 

S 2plXi 

T 2pl2i 

V 2pl3i 

V 2pl4i 

W 2pl5i 

Z 2pl6i 

T 2pl7i 

Z 2pl8i 

a 

2pl91 

b 2p20i 

c 

d 

• 

£ 

g 


h 

i 

j 

k 

1 

a 


n 

o 

P 

q 

r 

8 


t 

z 

u 

▼ 

w 

X 

y 


P 

1 2pl0m 

2 2plla 

3 2pl2m 

4 2pl3a 

5 2pl4x 

6 

2pl5a 

7 2pl6a 

8 2pl7a 

9 2pl8m 

A 2pl9a 

B 2p20x 

C 

2pl2o 

D 2pl3o 

E 2pl4o 

F 2pl5o 

G 2pl6o 

B 2pl7o 

X 

0 

a 

H 

3 

J 2pl9o 

K 2p20o 

X> 2pl4r 

X 2pl5r 

N 2pl6r 

o 

2pl7r 

P 2pl8r 

Q 2pl9r 

R 2p20r 

S 2pl6u 

T 2pl7n 

u 

2pl8u 

V 2pl9u 

W 2p20u 

Z 2pl8w 

T 2pl9w 

Z 2p20w 

a 

2p20y 

b 

c 

d 

• 

£ 

g 


h 

i 

1 

k 

1 

a 


n 

o 

P 

q 

r 

8 


t 

XL 

▼ 

w 

X 

y 


z 

1 2pl01 

2 2plll 

3 2pl21 

4 2pl31 

5 2pl41 

6 

2pl51 

7 2pl61 

8 2pl71 

9 2pl81 

A 2pl91 

B 2p201 

C 

2plln 

D 2pl2n 

E 2pl3n 

F 2pl4n 

G 2pl5n 

B 2pl6n 

X 

2pl7n 

J 2pl8n 

K 2pl9n 

L 2p20n 

X 2pl3q 

N 2pl4q 

o 

2pl5q 

P 2pl6q 

Q 2pl7q 

R 2pl8q 

S 2pl9q 

T 2p20q 

u 

2pl5t 

V 2pl6t 

W 2pl7t 

Z 2pl8t 

7 2pl9t 

Z 2p20t 

a 

2pl7v 

b 2pl8v 

c 2pl9v 

d 2p20v 

e 2pl9x 

£ 2p20x 

g 
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Not counting the extrapolation, up through n=20 there are 2382 even 
levels and 1962 odd levels, 182977 electric dipole lines between observed, 
Ritz, or polarization levels, and 286941 lines connecting levels predicted 
in the least squares fit. There are 36182 lines in series extrapolated from 
n«20 to ns99. There are 506100 electric dipole lines total. Forbidden 
lines were computed only for section X. There are 63223 magnetic dipole 
lines of which 16976 have good wavelengths. The are 398608 electric 
quadrupole lines of which 38124 have good wavelengths. Only 10 of the 
lines arise from metastable levels. 


If you do not wish to read about the calculation, the final results 
for the electric dipole lines and the 10 metastable forbidden lines are 
given in the files 


WL(NM) 
GF0600.POS 
GF0600EXTR. POS 
GF0 6 0 0 ALL • POS 
GF0600ALL. PRED 
GF0 6 0 0ALL • LINES 


WAVENO { CM- 1 ) 
GF0600 .WNPOS 
GF0 6 0 0EXTR . WNPOS 
GF0 6 0 0ALL • WNPOS 
GF0600ALL . WNPRED 
GF0 6 0 0 ALL . WNLINES 


NO LINES 


182987 

for 

36182 

for 

219169 

for 

286941 

for 

506110 

for 


"realistic" spectra 
line merging 
"realistic" spectra 
opacities 
opacities 


GF0 6 0 0ALL • LINES includes the other files. Predicted energy levels 
are listed as negative. In some files predicted lines are listed as 
negative. In those cases take the absolute value of wavelengths and 
energies. All wavelengths below 200 nm are vacuum and above are air. 
In the far infrared the wavelengths are essentially nonsense so use 
the wavenumbers. There is an annotated directory of the new files in 
0600 .README on the website. 


Transition integrals near nodes were found to be unreliable. The 
program that computes the transition array was modified to first read in 
all the Hartree-Fock transition integrals and then to read in 
scaled-Thomas-Fermi-Dirac integrals where they were thought to be more 
reliable. The current transition arrays were all Hartree-Fock. The 
program for computing the scaled-Thomas-Fermi-Dirac transition integrals 
has been rewritten to be at least an order of magnitude more accurate. 

Once the new calculations are systematically compared to other work, it 
will be possible to correct some of the transition integrals, if necessary, 
to improve the line list. 

Since I already had the NIST C I data compilation (Wiese, Fuhr, and 
Deters 1996) with 1300 lines in my computer, I made up tables comparing 
with my new line data for the allowed and forbidden lines, WFDK04.DAT, 
WFDK04BIGDIFF.DAT, FORBIDDENNISTK04.DAT The new data have intermediate 
coupling and configuration interactions which mean that I have to make a 
line by line check of transition integrals and eigenvectors before I can 
draw any conclusion about differences. I also compared with the lifetimes 
computed by Hibbert, Biemont, Godefroid, and Vaeck 1993), LIFEHBGVK04.DAT. 
My lifetimes are shorter by about 15% for lifetimes that "agree". Detailed 
comparisons will be made in the future. I will also compare computed 
spectra to the solar spectrum. 

1600. README 25 October 2004 1600 S I 

Added approximate Shore parameters for autoionization, gfl600.auto, 
and substituted gf 1600. auto into gfl600.pos to make gf 1600. sub. 
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20 October 2004 1600 S I 

The final results for tbs electric dipole, magnetic dipole, and 
electric quadrupole lines are in the files GFEMQ1600.POS for wavelengths 
in nm or in GFEMQ1600 .WHPOS for frequencies in wavenumbers in cm-1 for 
lines between observed levels. Only forbidden lines frost metastable 
levels are included. There are 2161 even and 2270 odd levels with 
configurations up to n * 15 or 16. There are 225605 electric dipole 
lines. Of these 24722 lines axe between known energy levels and have 
good wavelengths. The rest are predicted. 


From the hydrogen sequence data that have been published I have been 
able to interpolate os extrapolate the energy levels and to generate the 
lines up to n » 99 for el aments up through phosphorus. Isotopes are treated 
only through boron. The starting data for sulfur and higher are scant so 
that the interpolations and extrapolations may not be reliable. HIST should 
write a general program and compute all these data once and for all. 

I had an urgent request for hyperfine and isotopic splitting for Cu I so 
I incorporated the known splittings into my previously calculated Cu I line 
list and put the new data on the web. Fits to the solar spectrum look very 
good. 


I have copied ay extant test spectrum calculations to the web 
directory. The spectra were calculated either at resolving power of 500000 
or 2000000 in Doppler space so that a shift in point number is a Doppler 
shift. Dependi n g on effective temperature, the wavelengths range from 1 nm 
to 300$\mu$m. Same of the spectra have been broadened to lower resolution 
on the same point spacing and have been resampled to a smaller masher of 
points. Eventually the resolutions will include 100000, 50000, 30000, 
20000, 10000, 5000, 3000, 2000, 1000, 500, 300, 200, and 100. Oftentimes 
users can find what they need "off the shelf" on the web site. As I fill 
requests for new calculations, I add the new spectra to the website. This 
is the current list: alpha Can, Arcturus, A spectral types, Betelgeuse, 
Procyon, Sirius, Vega, GL411, Sun. 

Fiorella Castelli in Trieste and I used ay new distribution function 
program to compute new opacity tables and then new grids of model 
atmospheres with improved convection. The grids of models, fluxes, and 
colors have been put an the Kurucz web site. 

We have had complaints that we should be modern and use the Opacity 
Project and Iron project opacities. We (Castelli and Kurucz 2004) did an 
experiment in which she installed the Iran Project Fe I photoionization 
and autoionization opacity from Bautista (1997) into ay model atmosphere 
program ATLAS9 . The opacity looks huge with many autoionizing resonances 
in the continuum that required a great deal of supercomputer time to 
produce. In the Kurucz line lists the autoionizing resonances are treated 
< as line opacity on top of a smooth continuum opacity. I compute all the 

possible transitions between J levels while Bautista could treat only 
term-to-texm transitions. In some intervals I have much more opacity 
' because I have 39000 levels. For sense terms Bautista has much higher 

( smooth) continuum opacity and mine are wrong. The physical temperature 
structures of models computed with the two opacity treatments are the same, 
which implies that the total opacity, where it matters, is the same 
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averaged over small wavelength intervals for Bautista and Kurucz. In 
general the flux is the same as well except at one small interval where 
Bautista is higher and another where Kurucz is higher. Agreement does not 
imply that the models are correct. My old line lists are incomplete (that 
is why I am recomputing them) , so there must also be opacity missing from 
Bautista, and there are still systematic errors in the energy distributions 
predicted from the models. Comparisons with the observations (Bell, 
Balachandran, and Bautista 2001; Cowley and Castelli 2002) indicate that 
both computations appear to understimate the real Fe I absorption. 

And by implication all the iron group spectra. 

X have done minimal travel during this grant because it would slow 
down the work on the line data. 
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